Angular Dependence of Switching Properties in Single Fe Nanopillars 
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The continued increase in areal densities in magnetic recording makes it crucial to understand 
magnetization reversal in nanoparticles. We present finite-temperature micromagnetic simulations 
of hysteresis in Fe nanopillars with the long axis tilted at angles from degrees to 90 degrees to the 
applied sinusoidal field. The field period is 15 ns, and the particle size is 9 x 9 x 150 nm. The system 
is discretized into a rectangular pillar of 7 x 7 x 101 spins each with uniform magnetization. At 
low angles, reversal begins at the endcaps and proceeds toward the center of the particle. At ninety 
degrees, reversal proceeds along the entire length of the particle (save at the ends). The switching 
field was observed to increase over the entire range of angles, consistent with recent experimental 
observations. A second, lower-resolution micromagnetic simulation with 1 x 1 x 17 spins, does not 
agree with experiment, but shows behavior very similar to that of the Stoner-Wohlfarth model of 
coherent rotation. 
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Hysteresis in magnetic particles has been studied for 
nearly a century. However, the increasing areal densi- 
ties currently being pursued in magnetic recording tech- 
nologies are pushing bit sizes to length scales where 
the atomic nature of materials is important. In addi- 
tion, large data-transfer rates require switching times not 
much longer than the precession times of atomic mo- 
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FIG. 1: Hysteresis Loops for the 4949-spin model with mis- 
alignment angles of 0°, 45°, and 90°. The applied field is 
sinusoidal with a period of 15 ns, and the temperature is 
T — 100 K. One half of a loop was calculated and the data 
reflected to show a complete loop. 



ments. Taken together, further technological advances 
require a better understanding of hysteresis in nanoscale 
magnets at nanosecond time scales. Since these nano- 
magnets are often single-domain, the cooperative effects 
among the atomic spins determine the switching behav- 
ior. This cooperative behavior can be studied with sim- 
ulations of realistic model systems, and can be directly 
compared to high-resolution experimental resultsji*2i&4 
possibly for individual nanomagnets. 4 In the present pa- 
per we present such dynamic simulations of hysteresis 
in nanometer-sized Fe nanopillars exposed to a magnetic 
field that is aligned at an angle with the pillars' long axis. 

The general method in micromagnetics simulations is 
to divide a model magnetic system into coarse-grained re- 
gions, each with a uniform magnetization density M(fi) 
of a constant magnitude, Ms- Here ?i is the position 
of the ith spin in the system. The time evolution of 
each spin is given by the Landau-Lifshitz-Gilbert (LLG) 
equation 
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where H{?i) is the total local field at the ith position, 
7o is the gyromagnetic ratio (1.76 x 10 7 rad/Oe s), and 
a = 0.1 is a dimensionless damping parameter which 
determines the rate of energy dissipation. For the sign of 
the precession term, we follow the convention of Brown^ 
The total local field, H(fi), includes contributions 
from the applied field, dipole-dipole fields, and exchange 
interactions^ here the small crystalline anisotropy of Fc 



FIG. 2: The z-component of the magnetization at times t= 
4.48, 4.53, 4.58, 4.63, and 4.68 ns during the switching process 
of the 45° hysteresis loop in Fig. 1. 



is ignored. At nonzero temperatures, thermal fluctua- 
tions also contribute a term in the form of a stochas- 
tic thermal field which fluctuates independently for each 
spin. The thermal fluctuations are Gaussian with a 
mean of zero and covariance given by the fluctuation- 
dissipation theorem^ 

We have used two different models to simulate hys- 
teresis in single Fe nanopillars, as have been synthesized 
and studied by Wirth, et alm£££ The first model is a 
9 nmx9 nmxl50 nm pillar, discretized into a cubic lat- 
tice with 4949 sites (7x7x101). For such a large number 
of sites, direct calculation of the dipole-dipole interac- 
tions is impractical, and the Fast Multipole Method was 
employed, as described elsewhere^ The second model is 
a 5.2 nmx5.2 nmx88 nm Fe pillar discretized into a one- 
dimensional chain of 17 spins. This model assumes that 
deviations of the magnetization from uniformity over the 
cross-section of the pillar are negligible, as has been ob- 
served in Monte Carlo simulations i For both models, the 
material properties were chosen to correspond to bulk Fe, 
with saturation magnetization Mg=1700 emu/cm 3 and 
exchange length 2.6 nm. 

Figure 1 shows hysteresis loops at T = 100 K for the 
first model with the field misaligned at 0°, 45°, and 90° 
to the long axis of the pillar. The loops were calculated 
using a sinusoidal field with a period of 15 ns, which 
started at a maximum value of 10,000 Oe. In all the loops 
shown here, the reported magnetization is the component 



FIG. 3: The z-component of the magnetization at times t= 
2.25, 3.75, 5.25, 6.00, and 6.75 ns during the relaxation pro- 
cess of the 90° hysteresis loop in Fig. 1. 



along the long axis (z-axis) of the pillar. With the field 
and pillar aligned (0°), reversal initiates at the ends, as 
previously reported^ At 45° misalignment between pillar 
and field, the magnetization is initially pulled away from 
the long (easy) axis by the large magnetic field. As the 
field is swept toward zero, the magnetization relaxes until 
it essentially reaches saturation at zero applied field. 

Figure 2 shows the z-component of the magnetization 
at selected times during the reversal process for the 45° 
hysteresis loop in Fig. 1. As with 0°, regions of reversed 
magnetization nucleate at the ends of the pillar, and it is 
the growth of these reversed regions that leads to mag- 
netic reversal of the particle. It is important to note that 
the particles do not have a uniform magnetization, even 
though they are single-domain particles and would be in 
equilibrium under quiescent conditions. 

For 90° misalignment, the reversal mechanism is quite 
different. The hysteresis loop in Fig. 1 shows that the 
magnetization is essentially perpendicular to the easy di- 
rection until the field reaches a particular value. As the 
field is decreased further, the magnetization relaxes to- 
ward the easy axis. Since nothing breaks the up/down 
symmetry of the system when the applied field has no 
component along the easy axis, the relaxed magnetiza- 
tion can be directed toward either the positive or negative 
z-axis. Figure 3 shows the z-component of the magneti- 
zation for the 90° misalignment at selected times during 
the hysteresis loop. For this case, the relaxation occurs 
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FIG. 4: Hysteresis loops at 0°, 45°, and 90° misalignment 
between the field and easy axis for the 17-spin model. The 
applied field is sinusoidal with period 200 ns, and T = 10 K. 
Qualitatively, the loop shapes appear similar to those for the 
4949-spin model (Fig. 1), but the angular-dependence of the 
switching field for the two models is quite different. 
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FIG. 5: The angular dependence of the switching field _ff S w, 
shown the 4949-spin model (diamonds) and the 17-spin model 
(circles). The shape for the 17-spin model is qualitatively 
similar to what is expected from the Stoner-Wohlfarth model 
(dashed curve), while the 4949-spin model yields a depen- 
dence more consistent with recent experiments. 



along the entire length of the pillar, except at the ends. 
The pole-avoidance effect present at the ends, and in- 
volved in nucleation at smaller angles, retards relaxation 
along the easy axis. 

The hysteresis loops for the second, smaller model 
shown in Fig. 4 are qualitatively similar to those of Fig. 1. 
Again loops at 0°, 80°, and 90° misalignment are shown. 
There are important differences between the two models, 
however. First, without lateral resolution of the mag- 
netization across the cross-section, these pillars exhibit 
ringing due to the processional dynamics. Evidently, the 
precession of individual moments in the 4949-spin model 
does not lead to precession of the end-cap moment; pos- 
sibly the spin waves rapidly damp out the gyromagnetic 
motion. A second, and more prominent, difference be- 
tween the models is observed in the angular dependence 
of the switching field, H sw , shown in Fig. 5. Here H svf 
is defined as the applied field at which M z — (or in 
the case of 90° misalignment, where M z ^ 0). The sec- 
ond model (circles) shows a shape qualitatively similar 
to what is expected from Stoner-Wohlfarth (SW) theory, 
with a minimum H sw at 45°. The dashed curve is the SW 
theory with 7^ = 1600 Oe, much smaller than the 10 4 Oe 
assumed for ellipsoids. The first model (diamonds), on 
the other hand, has minimum H sw at 0°, and increases 
as the misalignment angle is increased. This behavior 
is consistent with recent experimental observations of Fe 
nanopillars i&l 

In summary, we have presented results from large- 
scale, finite-temperature micromagnetics simulations of 
hysteresis in models of nanometer-sized Fe pillars in mag- 
netic fields inclined with respect to the long axis of the 
pillars. In a model in which the magnetization is constant 
across the short dimension of the pillar, the angular de- 
pendence of the switching field is similar to that expected 
from Stoner-Wohlfarth theory. However, in a model with 
spatial resolution across the short dimension, the switch- 
ing field has a monotonically increasing angular behavior 
consistent with recent experimental results. 
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